Draft version 2013 February 4 

Preprint typeset using I^'T^]X style emulatcapj v. 5/2/11 



PRECISE HIGH-CADENCE TIME SERIES OBSERVATIONS OF FIVE VARIABLE YOUNG STARS IN 

AURIGA WITH MOST" 

Ann Marie Cody^'^, Jamie Tayar^-^, Lynne A. Hillenbrand^ Jaymie M. Matthews'*, and Thomas Kallinger'^ 

Draft version 2013 February 4 

ABSTRACT 

To explore young star variability on a large range of timescales, we have used the MOST satellite 
to obtain 24 days of continuous, sub-minute cadence, high-precision optical photometry on a field of 
classical and weak-lined T Tauri stars (TTS) in the Taurus- Auriga star formation complex. Obser- 
vations of AB Aurigae, SU Aurigae, V396 Aurigae, V397 Aurigae, and HD 31305 reveal brightness 
fluctuations at the 1-10% level on timescales of hours to weeks. We have further assessed the vari- 
ability properties with Fourier, wavelet, and autocorrelation techniques, identifying one significant 
period per star. We present spot models in an attempt to fit the periodicities, but find that we cannot 
fully account for the observed variability. Rather, all stars exhibit a mixture of periodic and aperiodic 
behavior, with the latter dominating stochastically on timescales less than several days. After removal 
of the main periodicity, periodograms for each light curve display power law trends consistent with 
those seen for other young accreting stars. Several of our targets exhibited unusual variability patterns 
not anticipated by prior studies, and we propose that this behavior originates with the circumstellar 
disks. The MO^T observations underscore the need for investigation of TTS light variations on a wide 
range of timescales in order to elucidate the physical processes responsible; we provide guidelines for 
future time series observations. 

Subject headings: circumstellar matter — open clusters and associations: individual (Taurus- Auriga) — 
stars:pre-main sequence — stars: variables: T Tauri, Herbig Ae/Be 



1. INTRODUCTION 

The widespread photometric activity seen in almost all 
young stellar objects is a defining characteristic of the 
class (Joy 1945, Herbig 1962). In particular, classical T 
Tauri stars (CTTS) display brightness variations froin 
the millimagnitude (iCodv fc Hillenbrai^ l20ia I2011D 
tomagnitude level ()Herbst et al.l 120021 iCarpenter et al.l 
I2001L |2002| ). on timescales from days to years. While 
some CTTS light curves appear to contain regular si- 
nusoidal patterns, the dominant form of the variations 
is aperiodic, with at times abrupt and unpredictable 
changes. The photometric variation of weak-lined T 
Tauri stars (WTTS; so called because their signatures 
of accretion are weak), on the other hand, appears to 
be dominated by spot rotation signatures, often with ad- 
ditional low level stochastic fiuctuations. These diverse 
properties have been variously attributed to rotational 
modulation of cool magnetic spots, enhanced chromo- 
spheric activity, hot spots from columns of magnetically 

|amc@ i pac. caltech.edu) 

^ California Instrfute of Technology, Department of Astro- 
physics, MC 249-17, Pasadena, CA 91125, USA 

Spitzer Science Center, California Institute of Technology, 
Pasadena, CA 91125, USA 

^ Department of Astronomy, The Ohio State University, 140 
W 18th Ave, Columbus, OH 43210, USA 

* Department of Physics & Astronomy, University of British 
Columbia, 6224 Agricultural Road, Vancouver, British Columbia 
V6T IZl Canada 

^ Institut fiir Astronomic, Universitat Wien, 
Tiirkenschanzstrasse 17, A-1180 Vienna, Austria 

Based on data from the MOST satellite, a Canadian Space 
Agency mission, jointly operated by Systems Canada Inc. 
(MSCI), formerly part of Dynacon, Inc., the University of 
Toronto Institute for Aerospace Studies, and the University 
of British Columbia with the assistance of the University of 
Vienna. 



channeled shocked infiowing gas, unsteady mass accre- 
tion, and occultation or shadowing by material in a sur- 
rounding circumstellar disk. Combinations of photomet- 
ric, spectroscopic and polarimetric techniques have illu- 
minated some of these possibilities, but a rigorous expla- 
nation of the erratic, often non-periodic nature of bright- 
ness fiuctuations awaits. 

Central to the problem of understanding photometric 
activity in T Tauri stars is the need for continuous mon- 
itoring on many different timescales. Photometry from 
ground-based facilities is limited in both precision and 
time coverage, consequently making efforts to model the 
light curves difficult, if not impossible. Routine gaps in 
observation are acceptable for the detection of periodic 
brightness variations, since these signals generally can be 
isolated in the frequency domain. But the fiickering por- 
tions of T Tauri lightcurves are not amenable to Fourier 
and similar analysis methods. 

To date only a few studies have monitored young stars 
at high precision and ca dence over long , conti nuous time 
baselines. Among these, lAlencar et al.l (|2010f) presented 
uninterrupted 23-day lightcurves of NGC 2264 members 
with 0.5-5 mmag precision from the CoRoT mission. 
They specifically identified objects with light curves re- 
sembling that of AA Tau, a CTTS with variability at- 
tributed to the repeated passage of warped disk mate- 
rial in front of the stellar photosphere every few days. 
The CoRo T dataset illustrated that AA Tau-like fading 
events are common in young stars with infrared excess 
and include not only dramatic brightness decreases but 
also erratic lower amplitude fiuctuations indicative of ad- 
ditional dynamics in the inner disk. The full complex- 
ity of this behavior was previously impossible to capture 
with groun d-base d time series. 

.Rucinski et al.l (pOlQi) used the Microvariability and 
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Oscillations of STars telescope (her eafter MOST; 
iWalker et al]|2003l iMatthews etalll2004[ ) to observe the 
Herbig Ae star HD 37806 for 21 days at ~3 mmag pre- 
cision, combining this dataset with nine seasons of ob- 
servations by the All Sky Automated Survey (ASAS). 
They characterized the light curve as weakly periodic 
on a timescale of ~1.5 days, and otherwise stochastic, 
with flares or accretion instabilities comprising ampli- 
tudes from 0.03% to 5% on timescales of minutes to 
ye ars. Another Herbig Ae star, HD 142666, was shown 
by IZwintz et al.l ()2009l ) to exhibit both 6 Scuti pulsa- 
tions and irregular UX Ori type variations attributed to 
the circumstellar disk. 

Additional high cadence investigation of young star 
photometric beh avior spanning multipl e we ek timescales 
was pre sented bv lRucinski et al.l (|2008D and lSiwak et al.l 
(|2011ai rbf). They acquired photometry at the 1-7 mmag 
precision level on TW Hya (nearly continuous observa- 
tions of 40 and 46 days) and five other T Tauri stars 
in the Taurus-Auriga and Lupus star forming regions 
(run durations of 12 or 21 days) from MOST. For the 
WTTS of their sample, most of the observed variability 
is well modeled by a collection of differentially rotating 
surface spots. The CTTS light curves present variabil- 
ity that is more challenging to interpret. In the case 
of RY Tau, there are two pronounced brightness dips of 
^^0.2 mag depth as well as lower amplitude, transient 
oscillatory behavior superimposed on a longer timescale 
trend. For TW Hya, on the other hand, flicker noise 
behavi or (power proportion al to inverse frequency) dom- 
inates (jRucinski et al.ll2008l ). although semi-periodic fea- 
tures were also observed to form and drif t to shorter 
timesc al es over the course o f a 40-day run (jSiwak et al.l 
l2011bD . ISiwak et all (|2011bD attribute this type of vari- 
ability to the magnetospheric accretion process, and in 
particular instabilities driving the flow of plasma from 
the inner disk. 

This small collection of space-based time series data 
has highlighted the complexity of young star variabil- 
ity and underscored the need for further datasets to de- 
termine the nature of the irregular and low-amplitude 
flux variations in these objects. The extent to which 
they are representative of young star light curves in gen- 
eral remains unclear. Currently the YSO VAR project 
(|Morales-Calder6n et al.]|2011l lRebul]||20Hl) is exploring 
the yield from multiwavelength precision photometry in 
several young clusters including the Orion Nebula Clus- 
ter and NGC 2264. 

We have taken advantage of the unique capabilities of 
the MOST mission to acquire a further high-precision, 
high-cadence, nearly continuous data stream over a 24- 
day observing period. Our aim was to monitor a hand- 
ful of erratically varying young stars to determine their 
short-timescale photometric patterns, decipher the mix 
of periodic and aperiodic phenomena in operation, and 
to quantify mathematically the aperiodic behavior. Pre- 
sented here is MOST data and analysis of the light curves 
of four T Tauri and Herbig Ae stars (SU Auriga, AB Au- 
riga, V396 Aur, and V397 Aur), and one new candidate 
early type Taurus member (HD 31305). 

The MOST observations are presented in Section 2 and 
the variability characterized in section 3. Section 4 con- 
tains discussion of the individual objects, and Section 5 
a general discussion of the implications of the photomet- 



ric timescales and amplitudes seen in our MOST 24-day 
observing sequence. 

2. OBSERVATIONS 

In operation since 2003, the 15 cm MOST telescope 
produces ultra-high precision differential photometry by 
virtue of its location 820 k m above Ea r th and specially 
designed imaging modes (jRowe et al.l 120061 IKuschni^ 
I2009f ). Since the satellite executes a polar orbit with 
period 101.413 minutes, it enables continuous viewing of 
targets for up to eight weeks in a zone (the "C VZ" ) from 
-18° to -1-36° declination. We selected a target field in 
the ~3 Myr Taurus- Auriga star-forming complex based 
on its location within the CVZ as well as the proximity of 
five suitably bright young targets: SU Auriga, AB Au- 
riga, V396 Aur, V397 Auriga, and HD 31305 (see Ta- 
ble 1). The former four are known young stars in the 
Taurus-Auriga association, while the latter is a newly 
suggested member based on this work. 

Observations took place over 24 days from 2009 De- 
cember 14 to 2010 January 07. All MOST images are 
acquired through a broadband filter with wavelength 
range 350-750 nm. Our three brightest targets (SU Aur, 
AB Aur, and HD 3130 5) w-ere m onitored in the direct 
imaging mode (Rowc e t al.l [20061 ) . which involves defo- 
cusing stars to a FWHM of 2-2.5 pixels. It can be ap- 
plied to up to 10 stars in the magnitude range U=6- 
11, and in our case produced photometry with point- 
to-point precision of 0.001-0.002 mag. For the two 
fainter ta rgets (V3 96 and V397 Aur), guide star imag- 
ing (i Walker et al.l |2005) offered the best performance, 
with prec isions of 0.01-0.02 mag. A third mode, Fabry 
imaging (jReegen et al.l l2006l ). produces photometry at 
the 10""* mag level but can handle only one target at a 
time and requires stars brighter than those in our sample. 

To mitigate pointing effects, MOST acquires and 
stacks many "subexposures" from 0.3-3 seconds. Total 
exposure times for the combined images was approxi- 
mately 30-40 seconds, with little dead time between; the 
resulting cadence was 43 seconds. Lightcurv es were gen- 
eratcd by the MOST ph otometric pipeline ()Rowe et al.l 
[2006, Hareter et al.ll2008D . which rejects cosmic rays and 
other artifacts, aligns the point spread functions (PSFs) 
for image stacking, and decorrelates background trends 
from target pixel fluxes. The latter stage was compli- 
cated by the high-amplitude variability present in most 
of our targets and hence required manual removal of clear 
outliers identified by eye. 

Following these reductions, there remained a resid- 
ual flux contribution from scattered Earthshine (e.g., 
iReeg en et all 120061 [Ro"we et al.ll2006D which introduced 
high levels of noise over a portion (~60%) of each satel- 
lite orbit. The data from these sections was unusable 
and hence not included in the final light curves. The re- 
maining flux points were modulated by a small (<0.5%), 
non-sinusoidal stray light contribution at the 14.2 cy- 
cles/day (165 microHz) orbital frequency. In addition 
to the once-per-orbit data gaps, there are also 15 larger 
gaps caused by observations of a source unrelated to our 
program. The resulting duty cycle is ~35%. 

Since our MOST data are not amenable to absolute cal- 
ibration, all photometry presented herein is relative. The 
resulting differential light curves all exhibit variability, 
such that the RMS values are not reflective of the under- 
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TABLE 1 
Targets and Basic Data 



Star 


Other identifiers 


B 


V 


SpT 


Noise RMS 


Empirical RMS 


AB Aur 


HD 31293, HBC 78, 2MASS J04554582+3033043 


7.1 


7.1 


A0/B9 


0.0017 


0.030 


HD 31305 


2MASS J04554822+3020165 


7.7 


7.6 


AO 


0.00092 


0.009 


SU Aur 


HD 282624, HBC 79, 2MASS J04555938+3034015 


10.2 


9.4 


G2 


0.0025 


0.110 


V396 Aur 


LkCa 19, HBC 426, TAP 56, 2MASS J04553695+3017553 


12.3 


11.2 


KO 


0.019 


0.050 


V397 Aur 


HBC 427, TAP 57NW, 2MASS J04560201+3021037 


12.9 


11.6 


K7I 


0.024 


0.058 



Note. — Basic data on the t arget stars, as repor ted by the SIMBAD database and measured from our time series data. Spectral types arc from 
SIMBAD, except for HD 31305 ||A rzner et al.ll2007h . Magnitudes are not simultaneous with our MOST time series observations. Noise RMS refers 
to the standard deviations of our MOST light curves after subtraction of smoothed median trend s, whereas the e mpirical RMS is the standard 
deviation inclusive of intrinsic variability. ^V397 Aur is a single-lined spectrosc opic binary system ( Walter et al.ll988f) : the spectral type corresponds 
to the primary. A direct imaging study determined AX— 0.87 for this system JKraus et al.|[201lD . 



lying white noise. We have estimated the point-to-point 
photometric precisions by subtracting out smoothed me- 
dian trends and fitting Gaussian profiles to the remaining 
noise distributions. The estimated precisions, along with 
the actual RMS values (before subtraction of the median 
trend), are listed in Table 1. 

3. VARIABILITY CHARACTERIZATION AND FEATURES 

The final MOST light curves are presented in Fig. 
1. Each of the five targets displayed variability well 
above the photometric noise level of a few mmag (Ta- 
ble 1). Brightness levels fluctuated by at least 0.03 mag 
(HD 31305) and up to 0.5 mag (SU Aur) over the course 
of the 24-day observations (see Table 1). In addition 
to the range of amplitudes, flux variation is present on 
a range of timescales. Periodic behavior predominates 
with characteristic timescales of a few days. Each star in 
our sample showed semi-periodic fluctuations. Further, a 
mixture of periodic and aperiodic components is evident 
in the light curves, in some cases with behavior that is 
not constant over the entirety of the light curve. For ex- 
ample, the last six days of the SU Aur light curve show 
a deep fading event; sir nilar episodes have b een reported 
in past literature (e.g.. iDeWarf et aT]|2003() . Several ob- 
jects display transient events, such as flaring. Fig. 2 il- 
lustrates this in two cases for the light curves of SU Aur 
and V396 Aur. 

In Section 3.1 we assess the periodic behavior using 
Fourier analysis. We then examine the evidence for corre- 
lated aperiodic behavior using wavelet analysis (Section 
3.2) and autocorrelation analysis (Section 3.3), which 
also affirms for the periodic objects the results of the 
Fourier analysis. 

3.1. Fourier analysis 

Because the light curve variations appeared to con- 
tain periodic components, we used Fourier analysis to 
identify frequencies and assess their significance. Among 
the unique features of this data were the long baseline 
and high cadence at which it was taken. The nearly 
uninterrupted sequence of MOST observations provides 
an opportunity to search for signals in the periodogram 
without interference from aliasing in the 1-10 day range 
typical of young st ar rotation and associa ted spot mod- 
ulation variability (|Irwin fc Bouvieii[2009l ). Over the 24- 
day duration, data points were generally spaced less than 
one minute apart. Additionally, the instrumental error 
on the brightness values was on the order of one millimag- 
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Fig. 1. — Differential light curves from MOST for the five targets, 
after removal of scattered Earthshine. The start of the observations 
(time=0) corresponds to Julian date 2455180.45. 
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Fig. 2. — Flares in two objects occurred during the MOST observations. 



TABLE 2 
Detected periodicities 



Star 


Periodogram 
period (d) 


Periodogram 
amplitude (mag) 


Wavelet 
period (d) 


Autocorrelation 
period (d) 


Previous values 
(d) 


AB Aur 
HD31305 

SU Aur 
V396 Aur 
V397 Aur 


6.551±0.015 
2.9379±0.0009 

2.661±0.001 
2.2393±0.0005 

9.54±0.02 
4.6022±0.0002 


0.0167±0.0003 
0.01056±0.00004 

0.0456±0.0002 

0.0652±0.0003 
0.04312±0.0004 

0.0477±0.0005 


5.97-7.55 
2.81-3.04 
2.39-2.64 
2.21-2.29 
9.36-9.50 
4.60-4.96 


o OO+0.03 
2.64±0.11 

9 r,.+0.05 

^■^^-0.04 
9 36+0.03 

4 44+0.19 


O.54I, 1.38^, 1.83^, 1.76^, 1.33* 

16^, 36, 1.7'^, 1.55*, 2. 73**, 2. 79,0.451° 
2.2311, 2.2412 
9.3213,9.3914^10.115 

4.7I6 



Note. — Periods, derived by identifying the largest periodogram peak and its 1-a uncertainties from Monte Carlo simulations (column 3), by 
identifying peaks in the wavelet function (column 4; range denotes the extent of the peak), and by selecting the first peak in the autocorrelation 
function and the timescales for which it drops by 1% from the peak value (column 5). The two rows under V 397 Aur correspond to its individual 
binary components. We list p re viously reported rotat io n periods for comp ari son. Reference s are a s follows: toorder et all l|2005l), ^atala et aT] 
11999). iTelleschi et al] II2007I). iPraderie et alj 11989). IPercv et al.l J2006I). 'IJohns fc Basril II1995I1 . 'IDeWarf eTall 120031) . IHerbst et al l 11987'), 
lUnruh et al.l 1120041). ^iFranciosini et al.l j2007[) . MGrankin et all l l2007l) .~^ 1Bouvier et al.l jl993bl) . ^ IZakirov et al.1 ngg^ .^ lNorton et al.1 {2007), 
i lBouvier et al.l J1995l ). ^Bouvier et al.l lll993aD 



nitude for the data on AB Aurigae, SU Aurigae and HD 
31305 and slightly higher for V396 and V397 Aur. The 
high quality of data allow us to search for periodic varia- 
tions ranging from multi-day rotation signatures down to 
minute-timescale brightness fluctuations, several mecha- 
nism s for which were pro posed recently (jKoldoba et al.l 
[200l lOrlando et al.ll20Tol) . 

Our selected tool for the analysis of periodic variabil- 
ity is the discrete Fourier transform (Deeming 197^ pe- 
riodogram. The MOST data sampling pattern includes 
a built-in periodicity as a result of the regular gap in 
observation (due to high levels of stray light; see Sec- 
tion 2) during the second half of each 101-minute orbit. 
Therefore, the periodogram of any non-constant bright- 
ness object contains an alias at the associated 14.2 cycles 
per day (d~^) frequency. We display an example of this 
phenomenon in Fig. |31 where a strong peak at /=14.2 
is visible in the periodogram of AB Aur. Because this 
peak repeats at 28.4 and other multiples of 14.2, we re- 
strict our initial period analysis to the frequency range 
from to 10 d~^ (i.e., periods greater than 2.4 hours). 
The Nyquist limit, or maximum independent frequency, 
is ~1000 d~^ (86 seco nd period). 

We used Period04 ()Lenz &: Bregeii I2005D to produce 
periodograms for each light curve. This program han- 
dles data with non-uniform time sampling and performs 



a multiperiodic least-squares fit for frequency, ampli- 
tude, and phase of a series of peaks selected in the pe- 
riodogram. We used Period04 in an iterative fashion: 
first searching for a single, most significant peak, then 
subtracting the corresponding best-fit sinusoid from the 
data, recomputing the periodogram, and searching for 
further peaks. We evaluated the significance of the iden- 
tified signals by requiring that their amplitude be at least 
a factor of 4.0 higher than the surrounding noise level in 
the periodog r am fo r 99.9% confidence, as suggested by 
iBreger et al.l (|1993D . In addition, a Monte Carlo simula- 
tion tool was used to add Gaussian noise to the best-fit 
combination of sinusoids, according to the noise level in 
the original light curve, and at the same time sampling. 
The periodogram was then regenerated, and the param- 
eters of the highest peaks recorded. We performed 500 
realizations of this process to determine the uncertainties 
on the periods and amplitudes quoted in Table 2. The 
full set of periodograms is displayed in Fig. 2] 

Using Period04, we discovered a single statistically sig- 
nificant period for each star (or two in the case of the 
V397 Aur system, wh ich is a spectroscopic binary that 
was recently resolved; iMathieu et al.l 119891 iKraus et al.l 
1201 If ), each of which is listed in Table 2. For V396 Aur, 
V397 Aur, HD 31305, the periodic pattern is relatively 
consistent over the 24-day duration of observations, and 
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Fig. 3. — A portion of the pcriodogram of AB Aur, demonstrat- 
ing the ahasing pattern at 14.2 that is associated with the 
MOST telescope's orbital 101-minute period. Power from intrinsic 
stellar variability in the 0—3 d~^ range leaks into regions of the 
periodogram with frequencies that are multiples of 14.2. 

the light curves encompass up to 10 cycles. We therefore 
ascribe the variability to stable hot or cool spots on these 
objects and tentatively associate the detected periodic- 
ities with their rotation rates. SU Atir also maintained 
stable periodic behavior consistent with spots, but only 
for the first ~18 days before a dramatic fading event set 
it. The case of AB Atir is even more; complex, since its 
light ctirve displays systematic trends on top of the sus- 
pected periodicity. The light curve of this star is too 
erratic to conclude that its brightness fluctuations arc 
reflective of its rotation period. 

We list the values of all detected periodicities in Table 
2. In several cases, they are consistent with one or more 
previously noted periods, suggesting that one or more 
spots are frequently present on the stellar surfaces and 
that the basic spot pattern is long-lived even if particular 
spots come and go over years. Nevertheless, the overall 
range of periods reported in the literature underlines how 
difficult it can be to infer accurate rotation rates from 
more sparsely sampled ground-based data. 

To assess even shorter timescales of variability, we ex- 
amined the full independent range of periodogram fre- 
quencies from to 1000 d~^ (or periods down to 1.5 
minutes). We performed this analysis on residual data 
after the sinusoidal signals detected with Period04 were 
subtracted out from the light curves in order to exclude 
power from the dominant periodicities. In addition, we 
masked out the MOST orbital aliasing features at multi- 
ples of 14.2 d-\ 

To model the frequency domain behavior, we fit 
a median trend to each periodogram and took its 
level at / ~800-1000 d^^ to represent the underlying 
white noise. These values ranged from 1.6x10^'^ mag 
(HD 31305) to 3.5x10--* mag (V397 Aur) and, for our 
dataset, are 2-5% of the photometric white noise in the 
time domain (i.e., the RMS light curve values listed in 
Table 1). We then fit a "1//" curve (i.e., power pro- 
portional to inverse frequency, or amplitude proportional 
to /-^/^) to the low-frequency end of the periodogram 
from 2-10 d"^ and noted the point at which this expo- 
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Fig. 4. — Periodograms, covering the 0—5 d~ frequency range. 
Vertical dashed lines mark rotation rates reported in the literature; 
in the case of V396 and V397 Aur, these closely match the position 
of our periodogram peaks. 

nential drop-ofi^ in variability amplitude with frequency 
reached the white noise limit. For AB Aur, SU Aur, and 

HD 31305, our targets with mmag-precision photometry, 
this occurred at ^^175 350 d^^. We detect no substantial 
deviations from flickering behavior (apart from the alias- 
ing power excess) on timescales down to several mimites. 
In addition, we can rule out high-frequency periodic sig- 
nals with periods down to 1 minute and amplitudes down 
to four times the periodogram noise limit: 0.4 mmag for 
AB and SU Aur, and 0.06 mmag for HD 31305. The re- 
sults of our Fourier analysis suggest that in general, ob- 
servations detect flicker noise in CTTS down to the white 
noise limit. Considering that ground-based photometry 
achieves a typical precision 5-10 times worse than our 
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MOST uncertainties, future observations would need tc 
sample light curves roughly once per hour to capture th( 
essence of variability in these stars. 

For the two lower precision datasets on V396 anc 
V397 Aur, amplitude levels in the periodogram foUowec 
a much shallower trend /~°'^^) after removal of th( 
main periodicities. We suspect that the light curves o 
these stars are fully explained by one or more spots alonj 
with low-level noise that is systematic but not character 
ized by flickering. We also fail to detect any significan 
high-frequency signals with amplitudes down to 1 mmaj 
in any of ou r targets. We concur from analysis of this 
sample with iGiinther et al.l ()201ClD . who found no evi 
dence for short-timescale periodicities due to oscillationi 
of an accretion shock. 

In summary, after removing the dominant periods du( 
to rotation and repeating the Fourier analysis, the result 
ing periodograms of the CTTS were all relatively feature 
less and consistent with: a "1//" flicker noise profile at 
low frequencies (<10 d~^), a white noise baseline at high 
frequencies, and a slight power excess due to aliasing at 
intermediate frequencies. Thus it appears that there are 
no characteristic timescales underlying the variability in 
these objects, other than the rotation period. A represen- 
tative example of the residual periodogram for AB Aur 
is shown in Fig. 5. 

3.2. Wavelet analysis 

While each of our target stars displayed a prominent 
signal in the periodogram, we found that most of the 
lightcurves are better described as a mixture of periodic 
and aperiodic phenomena. To determine whether some 
of the stochastic behavior could be attributed to the ap- 
pearance of coherent yet transient periodicities, we per- 
formed a wavelet transform using the WinWWZ pack- 
age available through the American Association of Vari- 
able Star Observer^]. This progr am applies a weighted 
wavelet Z-transform (jFosteii 119960 to achieve resolution 
in both frequency and time for unevenly spaced data. 
The wavelet employed is a simple sinusoid plus con- 
stant term, applied by sliding a window of predetermined 
width across the data. Datapoints close to the center of 
the window have the highest weight, whereas those near 
the edges are downweighted. WinWWZ is a useful tool 
for not only understanding how many periodic signals are 
present in the data at a given time, but also how their 
frequencies and amplitudes evolve. 

The basic output of the wavelet transform consists of 
signal power as a function of time and frequency. We 
illustrate this for our target stars with a series of con- 
tour plots, shown in Fig. 6. WinWWZ also provides 
one-dimensional plots of frequency versus time for the 
strongest signal. We detected one dominant period for 
each star (or two, in the case of binary system V397 Aur), 
affirming the results of Fourier analysis in Section 3.1. 
The derived values are presented in Table 2, with ranges 
indicating the minimum and maximum period values at- 
tained for this strongest signal over the course of the time 
series. 

The wavelet analysis suggests that there may be 
changes in period over the course of the MOST obser- 

WinWWZ was produced by G. Klingenberg and L. Henkel; it 
can be downloaded from http://www.aavso.org/winwwz 
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Fig. 5. — Logarithmic periodogram of AB Aur, after a 6.5 day 
period is removed from the light curve. We have fitted a "1//" 
profile (i.e., amplitude proportional to f~^^^) to the low-frequency 
regime (dashed diagonal line), and a uniform noise level at high 
frequencies (dashed constant line) . The series of peaks starting at 
log(frequency)=1.15 are aliases due to the MOST satellite orbital 
period; the first seven are marked with vertical ticks. 

vations. The abrupt disappearance of the periodic sig- 
nal in the light curve and wavelet transform of SU Aur 
around day 18 is the most obvious instance of a change 
in periodic behavior. Since we do not know a priori the 
shape of the subsequent brightness dip, it is difficult to 
ascertain whether the periodicity vanishes completely or 
is simply washed out by the strongly systematic light 
curve trend. If the former, then the disappearance of 
the signal indicates that the stellar photosphere becomes 
obscured by material in the line of sight. Other more 
subtle variations in period include a 1.5 day evolution in 
the strongest period of AB Aur, and a change of over 0.2 
days in the period in the HD 31305. These effects may 
reflect actual changes in the variability timescale, or in- 
stead be attributed to non-sinusoidal components in the 
variability. 

As the target with the most stable frequency behavior, 
V397 Aur shows two strong and sustained periodic sig- 
nals in the wavelet analysis, consistent with its identifica- 
tion as a binary system. For the other objects, the ranges 
in period measured for the dominant signal are of order 
100 times larger than the l-cr uncertainties returned by 
the periodograms. Since the periodogram represents the 
average behavior of periodicities over the 24-day run, the 
wavelet result again suggests that there is subtle evolu- 
tion in the timescales of periodic variability during our 
observations, or that aperiodic variability is interfering 
with out ability to measure a stable period. 

3.3. Autocorrelation analysis 

Fourier analysis and wavelet analysis provided an in- 
dication of the periodic vari ations present in the data. 
Autocorrelation (|Box fc Jenk ins 1976) is a further tech- 
nique that enables confirmation of these results as well 
as additional searches for variability that is not strictly 
periodic, not sinusoidally shaped, or not present for the 
entirety of the data set. It provides an assessment of how 
consistent variability patterns are on different timescales. 
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Fig. 6. — Wavelets for all targets, zoomed in on the frequencies of significant signals. The intensity scaling is arbitrary but the same for 
all objects. The Fourier periods are recovered, but some additional structure is also apparent. 
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Perfect correlation results in an autocorrelation value of 
1, while completely uncorrected data returns a value of 
0, and anticorrelated data corresponds to negative val- 
ues. This technique is particularly useful in that it can 
reveal periodic variations even if they are not persistent 
throughout the entirety of the data set. It serves as a 
check on the results obtained by Fourier analysis and 
wavelet analysis. Since the MOST data were not evenly 
spaced throughout the observation period, we first re- 
sampled the light curves using linear interpolation. We 
then applied the unbiased autocorrelation formula: 

^ Jmax 

\y /Jmax ^.^j^ 

where N is the total number of data points, At the 
timeshift, y the light curve values, and jmax = N — t/At. 

We noted a single prominent autocorrelation peak 
for each star, consistent with the periods derived from 
Fourier analysis. We have derived periods by measuring 
the timeshift of the first maximum and noting the sur- 
rounding values for which the autocorrelation function 
drops by 1%; these ranges are listed in Table 2. Since 
the light curve of SU Aur ceases to display oscillations 
during the last ~6 days of observation, we have computed 
two versions of the autocorrelation function for this star. 
The first is performed on the entire dataset, while the 
second includes only the first 18 days. 

In addition to providing confirmation of the previously 
measured periods, autocorrelation analysis enables inves- 
tigation of any further underlying variability. To accom- 
plish this, we first used Period04 to subtract out vari- 
ability at the dominant period. We then recomputed the 
autocorrelation function on this residual. Both versions 
of the autocorrelation are shown in Fig. 7. The residuals 
contain low-level undulating features, but further sub- 
traction of periodogram-detected signals from the light 
curve did not result in their disappearance. 

While no significant new variability signals were evi- 
dent in the residual autocorrelations, the analysis does 
provide an indication of the characteristic timescales for 
aperiodic physical processes in our target stars. We as- 
sociate an autocorrelation timescale by noting the value 
for which the autocorrelation drops to 0.5. The flux vari- 
ations of HD 31305 are correlated on timescales of ^0.2 
day once the periodic signals are removed, while SU Auri- 
gae is correlated on timescales of ~1.1 days. In contrast, 
AB Aurigae displays coherence in its variations out to 6.2 
days. Long timescale correlation after the periodic sig- 
nal has been removed implies that for these three sources 
there are additional physical processes other than peri- 
odic spot modulation that contribute to the variability. 
The residual light curves of the two weak-lined T Tauri 
stars V396 and V397 Aur are correlated on much shorter 
timescales of <0.1 days. We suspect that spots explain 
most of the variability in these systems. 

4. STARSPOT MODELING 

The periodogram, wavelet, and autocorrelation analy- 
sis methods yielded one significant periodicity per star. 
In most cases, the derived timescale was very close to 
a rotation period previously reported in the literature. 
The long-term stability of these periodicities, as well as 
the deviations in light curve shape on shorter time scales. 



suggests that the underlying mechanism is flux modula- 
tion by cool or hot starspots, as opposed to stellar pul- 
sation or an eclipsing companion. 

Representation of a light curve with spots requires a 
number of parameters and often admits degenerate so- 
lutions. A non-parametric model involving only total 
spot coverage and temperature contrast relative to the 
photosphere is one way to avoid these challenges (e.g., 
iGrankin et al.l [20081 ). However, the dense sampling and 
long duration of our time series observations enables 
more detailed exploration of the presumed spot prop- 
erties than is typically possible with ground-based data. 
We exploit the high quality of our datas et by applying 
StarSpotz ()Croll et al.ll200a iWalker et al.ll2007[ ). an ana- 
lytic starspot modeling program designed specifically for 
the analysis of MOST data. 

To reduce the number of free parameters required for 
the spot modeling process, we fixed the stellar inclina- 
tion, i, as well as the rotation period. We set the latter 
according to our results from Fourier analysis (Section 
3.1). The value of i was taken from the literature where 
available. In addition, we fixed the limb darkening coef- 
ficient via the l/-band prescription of Claret & Blocme^ 
(|2011D . but tested several alternate values to determine 
the sensitivity of our results to this parameter. 

We provided starting guesses for all other free param- 
eters, including unspotted stellar brightness, spot dark- 
ness (ratio of flux compared to that of the surrounding 
photosphere), spot latitude and longitude, as well as spot 
angular size. To avoid unphysical values of spot darkness, 
we produced models for discrete steps between 0.0 (com- 
pletely black spot) and 1.0 (same brightness as photo- 
sphere). The total number of spots, as well as the pres- 
ence of differential rotation, were specified in advance, 
and the free parameters were then varied separately for 
each spot to create a suite of analyt ical models fo l lowin g 
the prescriptions of iBuddin^ ()1977D and iDorrenI ()1987D . 
Starspotz evaluates the goodness of model fits to the light 
curve via a non-linear least squares algorithm, returning 
a local minimum value. 

We assessed the improvement in fit brought about by 
increasing the number of spots, as well as whether an 
assumption of differential rotation produced a better fit. 
In general, neither of these additions resulted in signifi- 
cant improvement and there remained deviations in the 
light curves that could not be explained by spots. Since 
these additional features are lower in amplitude than the 
periodic variability, we argue that spot models can help 
nevertheless to inform our picture of stellar surface prop- 
erties. In the following sections, we detail the the best-fit 
spot configurations and evaluate their consistency with 
the light curves of V396 Aur, SU Aur, HD 31305, and 
AB Aur. We did not produce a spot model for V397 Aur 
since it is a binary system and both members appear to 
contribute to flux variations in the light curve. 

4.1. V396 Aurigae 

Several rotation periods for the star V396 Aur were 
previously report ed in the literature. T he 2.228 days 
by Granki n et al.l ((2008) and 2.24 days bv lBouvier et al.l 
(fl993bl agree well with our values derived from Fourier, 
wavelet, and autocorrelation analysis. The former study, 
known as the Maidanak survey, characterized variabil- 
ity in this object as having a correlation between the 
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Fig. 7. — Autocorrelations (solid curves) and residuals (dashed) after period removal. We have used the residuals to estimate a char- 
acteristic timescale for variability, by noting where the autocorrelations drop to 0.5. The second panel for SU Aur is the autocorrelation 
performed on only the first 18 days of the time series, before the prominent fading event begins. 



magnitudes of maximum and minimum light during dif- 
ferent observing seasons. An explanation for this be- 
havior could be long term evolution of the total spot 
area, which they represent by a collection of cool spots 
on the stellar surface. Despite the history of observation, 
the typical distribution of the starspots causing variabil- 
ity in V396 A u r rem ains ill determ ined. In cont r ast to 
iGrankin et all (|2008| )'s conclusions. iHuerta et all ()2008f ) 
proposed a single large spot near the star's pole, based on 
2.18-day periodic radial velocity variations of ~1 km s^^ 
amplitude. 

We set out to model the starspot configuration of 
V396 Aur during our MOST observations by fixing 
known stellar parameters and generating a series of one 
to three dark spots. The vsmi value of V396 Aur is 



18.6±1.9 km s-i (jHartmann et al.lll987l ). yielding stel- 
lar inclinations from ^^25-35° for the range of range 
of radii expected f or a young KO star {R ~1.5 Rq; 
iBouvier et al.l I1995D . To generate Starspotz models, 
we adopted the two r nost extreme values reported in 
the literature : i=3 3° (|Bouvier et al.l Il995| ) and i=25° 
(jHuerta et al.l 120081 ) . along with a limb darkening coef- 
ficient u = 0.66 appropriate f or the KO spectral type o f 
V396 Aur (Toff ~5250 K) from lClaret fc BloemenI (|20T1 . 
Combined with the adopted vsmi and stellar radius, our 
observed photometric period of 2.24 d is consistent with 
values of i in the range 33°±4. 

We find that the MOST observations of V396 Auri- 
gae are best modeled by a single large spot at 60-80° 
latitude, with the lower latitudes favored for the stellar 
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Fig. 8. — Left: Lightcurve of V396 Aur, binned on the MOST orbital timescale (points with uncertainties), compared to one of the best 
fitting spot models (solid curve). Right: Unbinned light curve residuals after model subtraction. 



inclination i=25°. The best fitting solution is dependent 
on the spot darkness parameter, with darker spots corre- 
sponding to smaller covering fractions. We found a min- 
imum spot diameter of 23°, or ~4% of the stellar surface 
area). We find acceptable values for the spot darkness 
between 0.0 and 0.5; spots with lower contrast would oc- 
cupy an unphysically large portion (>30%) of the visible 
photosphere. The solution is also degenerate with the 
unspotted stellar brightness, again leading to very large 
spots for unspotted brightnesses much greater than the 
observed flux level. Incorporation of two or more spots 
to the model did not result in a significantly better fit 
to the light curve. Adjustments to the limb darkening 
coefficient also failed to bring about improvement. 

In Fig. [8] we overplot one of the best fitting models 
on the light curve of V396 Aur, with data binned on 
the 101-minute MOST satellite orbital timescale, along- 
side the residuals of the unbinned light curve after model 
subtraction. There are a number of outliers from the 
model, one of which occurs at i=3 days into the run. 
Closer inspection of the full cadence light curve reveals 
that this is a flare event (see Fig. 2). The rise and de- 
cay times are 0.5 and 1.2 hours, respectively, and the 
amplitude is '^0.1-0.15 mag. Other outliers lie mainly 
at the peaks and troughs of the light curve, suggesting 
that spot evolution is taking place, or that there is an 
additional source of low-level variability that is most ap- 
parent at these phases. Nevertheless, there do not appear 
to be any systematic trends or periodicities in the residu- 
als. The minimum reduced value was 36.1, confirming 
that a spot model does not provide a full explanation for 
the variability in V396 Aur. 

4.2. SU Aurigae 

Previous observation of the G2 star SU Aurigae has in- 
dicated periodic variability on a variety of derived time 
scales (Table 1). A number of observers have also noted 
its erratic dimming b y ^0.5 mag in the o ptical over time- 
frames of 3-5 days (jUnruh et all I2004D . This has led 
to its classification as a UX Orionis star ("UXOR"), 
although the change in magnitude is not as large as 



that of prototypical UXORs (^1 mag in V\ Herbst & 
Schevchenko 1999) or slightly lower mass "Type III" T 
Tauri stars (<1 mag; Herbst et al 1994). 

Our observations of SU Aurigae also indicate a picture 
of variability involving a mixture of behavior. We observe 
both a periodic component during the first 20 days of 
observation and an 0.4 mag dimming event that begins 
around i=19 d and persists through the end of the time 
series. In addition, a flare occurs at i 10.2 d, with rise 
time '^1 hour, decay time ~2.4 hours, amplitude 0.03- 
0.04 mag (see Fig. 2). At i 1.7 d, a portion of the light 
curve appears offset and brighter by ^^0.04 mag. As far 
as we can tell, this is not an instrumental error, and the 
star was inexplicably brighter during this particular peak 
of its oscillatory pattern. 

To model the periodic variation in SU Aur, we have 
considered only the first 20 days of the time series, 
since the last five days are dominated by systematic fad- 
ing. Our chosen limb darkening c oefficient is u = 0.63 
based on iClaret fc Bloemeiil (|2()11D. and appropria te for 
a young G2 star (Toff - 5550 K: IDeWarf et"an[200l . We 
chose a stellar inclination by assuming that it is equal to 
the inclination measured interferometrically for the disk: 
i = 62+g° (| Akeson et all |2003) . The spectroscopic ro- 
tation v elocity has been consistently measured at ^66 
km s~^ (iBouvier et al |[l986llHartmann fc Stauffeiill989l 
iGlebocki fc Gnacihskil l20()5l ). We ad opt the value of 
wsim=66.2±4.6 km s""'^ reported by iHartmann et al.l 
(1986). Combining these measurements with values of 
the steflar ra dius reported i i i the literature {R = 2.75 ± 
0.25-R,n from iDeWarfet all (|200l and R ~3.6 Rq from 
iPetrov et all (|1996( )). we find a range of expected rota- 
tion periods from 1.9-2.4 days. Including propagation of 
uncertainties in inclination and rotational velocity, this 
range expands to 1.7-2.6 days. The inferred periods are 
marginally inconsistent with the 2.66d timescale that we 
have measured photometrically. 

Assuming that the variability is caused by features on 
the stellar surface, we find that the periodic component of 
the light curve could be accounted for by either a dark or 
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Fig. 9. — Left: Lightcurve of SU Aur, binned on the MOST orbital timescale (points with uncertainties), compared to one of the best 
fitting spot models (solid curve). Right: Unbinned light curve residuals after model subtraction. 



bright spot at high latitude (80-85°). For the dark spot, 
we determine plausible ratios of spot to photospheric flux 
between 0.0 and 0.3, and a spot size of ~28-40° (~6-12% 
of the total stellar surface area). A bright spot, on the 
other hand, could be a factor of 2.0 or more brighter 
than the surrounding photosphere and up to 30-40° in 
size (7-12%). 

As with V396 Aur, neither the addition of more spots 
nor adjustments in the limb darkening parameter offers 
substantial improvements in the model fit. We note that 
with the single band MOST observations, we are unable 
to distinguish a hot spot generated by an accretion shock 
from a cool spot generated by the stellar magnetic field. 
However, the unstable amplitude of variability from one 
cycle to the next in our light curve suggests that ac- 
cretion may be a better explanation. Since the model 
light curves appear nearly identical, we overplot only the 
dark spot model on the orbit-binned light curve in Fig. 9, 
alongside the residuals after model subtraction. The full 
light curve of SU Aur shown in Fig. 1 displays inconsis- 
tent amplitudes from one cycle to the next. Examination 
of the residuals in Fig. 9 reveals that there are both pos- 
itive and negative deviations, with no systematic trend. 
The minimum reduced value is over 4000 for all mod- 
els, indicating a poor fit. We therefore suspect that there 
is an additional source of variability present, or that the 
periodic mechanism involves dynamic evolution on '^l 
day timescales. In Section 5.3 we propose that orbit- 
ing material at the inner disk edge of SU Aur offers an 
additional explanation for the observed periodicity. 

4.3. HD 31305 

HD 31305 is not a known member of the Taurus- Auriga 
star-forming region, but we include it in our presentation 
since it is also a bright star in our MOST field, and it 
exhibited surprising periodic variability. Full interpreta- 
tion of the variability observed with MOST hinges on the 
status of this AO target as a potential young star, a ques- 
tion which we defer to Section 5.2. A series of pronounced 
cyclic flux changes suggest that either an inhomogeneous 
stellar surface or orbiting circumstellar disk features play 



a role in the variability. Stellar pulsation is another peri- 
odic phenomenon present in early-type stars, but obvious 
inter-cycle changes in the light curve do not support this 
alternative explanation. With a period near 2.8 days and 
amplitude of '^0.01 mag, the variability pattern displays 
significant deviation from sinusoidal behavior at the 0.5% 
level, none of which corresponds to significant periodici- 
ties in the periodogram. 

We have attempted to model the overall behavior of 
the light curve with a series of dark or bright starspots. 
These could represent either regions of cooler surface 
temperature, patchy photospheric chemical composi- 
tion s uch as observed in s ome peculiar typ e A stars 
(e.g., IShulvak et aP [20lol IKochukhovl mm . or hot 
spots related to low-level accretion. The inclination of 
HD 31305 is unknown, so we estimate it by consider- 
ing the stellar parameters in concert with the vaim. value 
of 50±25 km s""'^ measured by Mooley et al. (2013, in 
preparation). Of note, this vsmi puts HD 31305 at the 
low end of the distributi on of rotation rates for A stars 
fsee lZorec fc RoveHlMl . , Fig. 7). Using U, B, V, J, H, 
and K magnitudes from the Simbad database along with 
a spectral type of AO, we derived a bolometric luminos- 
ity based on SED fitting. Using only the optical magni- 
tudes, we derive a luminosity of log(L/L0)=1.38. Since 
there is an infrared excess associated with this object, 
non-photopheric emission likely contaminates the SED 
at J band and beyond; incorporating the magnitudes 
here results in a larger luminosity of log(L/LQ)=1.5. 
This range of luminosities, along with an AO effective 
temperature of '^9850 K results in a stellar radius of 
1.85±0.09 Rq. If the observed periodic variability is due 
to surface starspots, then the combination of the 2.97 d 
period and estimated radius require a t;sinz of at most 
31 km s~^. For a minimum vsini of 20 km s~^, we de- 
rive i=39°. To perform the spot modeling, we therefore 
adopted representative values of 40° and 70° . 

From the configurations generated with Starspotz, we 
find that the light curve can be explained by either a 
single high-latitude (70-90°) spot subtending 1-20% of 
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Fig. 10. — Left: Lightcurve of HD 31305, binned on the MOST orbital timescale (points with uncertainties), compared to one of the best 
fitting spot models (solid curve). Right: Unbinned light curve residuals after model subtraction. 

curve, SO that models could be fit to the shorter timescale 
(<10 d) variability in AB Aur. We adopted two values 
for i: 22° and 3 5°, matching the disk inclinatio n val- 
ues presented by I C order et al.l (2005) and P errin et alj 
12009). Assuming one mid-latitude (20-60°) spot sub- 
tending 10-15° in diameter, we could reproduce some of 
the major trends of the data (see Fig. 11), and with a 
large number of small spots, we also reproduced some of 
the smaller fluctuations. The results were relatively in- 
sensitive to the limb darkening parameter and the spot 
darkness. As suspected, none of the spot coverage sce- 
narios provides a particularly satisfactory model, with 
reduced values in excess of 3000. However, we can- 
not rule out that a spotted photosphere contributes to 
a portion of AB Aurigae's variability. As discussed for 
HD 31305, the light curve behavior observed here is not 
readily understood in an AO star. 

5. SUMMARY AND DISCUSSION 

We have presented high precision, high cadence light 
curves for four known Taurus- Auriga members as well 
as one newly suspected member, HD 31305. The stellar 
masses for these stars range from ~0.5-2.5 Mq. While 
our sample is not large, the high quality of the data 
has enabled us to probe young star variability down to 
sub-millimagnitude levels and search for low-amplitude 
and/or short-timescale phenomena. 

5.1. Discussion of lightcurve analysis 

Our investigation has resulted in the measurement of 
photometric periods for these objects, all of which dis- 
play a single dominant periodicity. For cases in which 
the periodic behavior can be explained by starspots, we 
have modeled the possible configurations and found that 
one mid- to high-latitude spot provides a sufficient, but 
far from perfect, explanation for the periodic variabil- 
ity. This pic t ure is largely consistent with the results of 
iSiwak et al.l (|2011a[ ). who carried out spot modeling of 
several variable young G and K stars. Among the pe- 
riodic stars of their sample, they generally find a high 
latitude or polar spot (often large), and sometimes an 



the stellar surface or two spots at a large range of lati- 
tudes, with diameters from 3-20° (up to 3% of the area). 
The best fitting two-spot solutions depend highly on the 
adopted unspotted stellar brightness, as well as the spot 
darkness, and do not favor any particular latitude. We 
attempted to add three spots, but this did not improve 
the fit. 

We have overlaid one of the best fitting models for 
HD 31305 on its light curve in Fig. 10. Deviations at the 
0.5% level, as well as a reduced value larger than 1100, 
confirm that there remains significant aperiodic variation 
in the data. 

4.4. AB Aurigae 

AB Aurigae is a Herbig Ae star with spectral type 
AO. A large body of observations has indicated that it 
is encircled by an optically t hick, gas dominated disk 
from 0.24 AU to 300 AU (e.g.,[Ta nnirkulam et"an [200l 
IQppenheimer et al]|2008f) . This disk is viewed almost 
face-on, at an in c hnation angle between 12 and 35° 
(iPietu et all 120051 ^agawa et al.l 12004 iCorder et akl 
l2005llMarifias et alll2006l) ! 

Previous photometric studies of AB Aur uncovered pe- 
riodic variability on short timescales from 0.5-1.8 d (Ta- 
ble 2). Our analysis found a much longer period (6.59 d), 
with substantial aperiodic light curve components on 
both longer and shorter timescales. Evidence that this 
is not the star's rotation period comes from the reported 
spectroscopic rotation velocity of wsim=80±5 km s~^ 
(|Bohm fc Catalal 119931 ) and larger, wit h a radius of 
2.0-2 5 assuming L=38-4 7 from [PeWarf et all 
(|2003l ) and iTannirkulam et all (HOM)- These values are 
inconsistent with such a long period; even if i were as 
large as 90°, the photometric period would have to be 
less than 1.6 days. 

Consequently, the variability observed in AB Aur is not 
particularly amenable to spot model fitting. For com- 
pleteness, and in case any of the reported parameters 
are erroneous, we have proceeded nevertheless to model 
its light curve using Starspotz. To begin, we subtracted 
out two linear trends at days 0-5 and 5-25 in the light 
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Fig. 11. — Left: Lightcurve of AB Aur, binned on the MOST orbital timescale (points with uncertainties), compared to one of the best 
fitting spot models (solid curve). Right: Unbinned light curve residuals after model subtraction. Linear trends have been removed from 
the left light curve, but not the right. 



additional smaller spot closer to the equator. 

The main difference between our results and theirs is 
that we are unable to account for the full light curve be- 
havior in our periodic objects using spot models alone. 
Additional aperiodic variability appears to be prevalent 
in the higher mass stars of our sample (AB Aur, SU Aur, 
and HD 31305) at the 0.1-1% level, based on the flicker 
noise profiles in their periodograms and the substantial 
residuals after subtraction of best-fit spot models. The 
variability patterns among these three objects is never- 
theless distinct, with differing characteristic time scales 
found by our autocorrelation analysis (Section 3.3). In- 
vestigation of a larger sample of young variable stars 
may shed light on how these differences arise; we sus- 
pect that the inclination of surrounding disks and ac- 
cretion columns plays a substantial role in the prop- 
erties of observed variability. In t he case of SU Aur 
the dust is viewed ne arly edge on (jAkeson et al.l [20021 
[Kurosawa et al.l 1200 5*). while for AB Aur, it is close to 
face on (jMarifias e t al. 2006). An other nearly face-on 
object, TW Hya, was shown to by ISiwak et al] (|2011bf ) 
to exhibit intermittent and variable-period oscillations 
mixed with more stochas tic variability. MHD si mula- 
tions of disk accretion by iRomanova et al.l (|2008l ) have 
shown that there is a regime of suitably high accretion 
rates for which the flow of material becomes unstable, 
and the resulting light curve is stochastic. 

The accretion rate of SU Aur has be en measured 
at a moderate 0.5-0.6 x IO'^Mq yr"! (jCalvet et al.l 
120041 ). Accretion rates reported for AB Aur, on the 
other hand, vary widely, from re latively high values 
of 1.4 and 3xlO~^M yi"^ ( Garcia Lopez et all [20061 
IBrittain et all I2007D to 7 x 10' ^ Mp^ yr"^ and lower 
based on more indir ect methods (jBohm fc Catalal [19931 
iTelleschi et"alll2007[) . It is therefore difficult to surmise 
whether how the variability observed in our target stars 
is connected with accretion properties. 

Finally, we note that this work has resulted in the mea- 
surement of photometric rotation periods for both mem- 
bers of the binary V397 Aur. One period is close to but 



not exactly twice the other, suggesting that a tidal reso- 
nance may be operating in this system. The detection of 
multiple periods in a combined light curve also highlights 
the possibility that more close young binaries could be 
detected in this way without the need for extensive ra- 
dial velocity data. For binaries in which the two rotation 
periods are similar, care must also be paid to distinguish 
the scenario from differential rotation. 

5.2. HD 31305: An unusual variable young A star? 

As discussed in Section 4.3, the MOST light curve of 
HD 31305 resembles that of an active young star. Only 
recent analysis (see below) has hinted that it may be a 
member of the few-Myr-old Taurus- Auriga star-forming 
complex. All prior literature has assumed it to be a 
background field star, and in some cases, used i t as a 
photometric comparison star (e.g., iHerbst et al.l I1987L 
[PeWarf et al...2003,') ! 

A clue to its potential youth status is the infrared 
excess detected by IRAS (identifier 04526-^3015), im- 
plying that HD 31305 is encircled by dust at ^390 K 
(jOudmaijer et aL 1992^ Recent observations from the 
WISE mission (RebuU et al. 2011) corroborate the pres- 
ence of infrared excess, although the star was listed 
as a background object due to the absence of Taurus- 
Auriga membership informatio n. It star has also been 
reported to fiare in the X-ray (lArzner et al.ll2007f) with 
a peak temperature of 8.6 keV (Franciosi ni et alT ^OO?!. 
These properties are largely inconsistent with those of A- 
type main seq uence stars, which are weak X-ray emitters 
(|Linskvl[2003D . Additional evidence from Mooley et al. 
(2013, in preparation) bolsters the idea of a young age for 
HD 31305. They report a spectroscopic distance (174±11 
pc) and proper motion consistent with Taurus member- 
ship at the 80% confidence level. Their assembled SED 
from optical through mid-infrared wavelengths displays 
indications of a 10 /lui silicate feature, and a model fit 
implies a temperature of ^^350 K at the inner edge of 
the dust region, consistent with the value suggested by 
.Oudniaiier et al., (,1992i) . However, their low- resolution 
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spectrum lacks emission li nes, suggesting that t h e disk 
is not actively accreting. iMiroshnichenko et alj (|1997[ ) 
pointed out that infrared excesses in HD 31305 and other 
A-type stars may be consistent with circumstellar dust 
shells. We alternatively suggest that it is weak enough 
to implicate a debris disk, since the ratio of infrared to 
stellar luminosity in the SED assembled by Mooley is 
~1:200. This scenario would not be inconsistent with an 
A star that is significantly older than ^10 Myr, rather 
than the 3 Myr age of Taurus. Therefore, we cannot rule 
out that HD 31305 is a main sequence star, although 
its light curve behavior remains atypical and difficult to 
explain. 

The variability itself has features in common with 
the other stars in our sample such as AB Aur, specifi- 
cally a light curve including both periodic and aperiodic 
components. To our knowledge, such a mixture of be- 
havior lacks a mechanism among main sequence stars. 
The most common explanation for non-eclipsing periodic 
variability- magnetic surface spots- is difficult to invoke 
for A stars since they are not believed to host substan- 
tial magnetic fields. The convective/radiative boundary 
vanishes for effective temperatures above ^6000 K, im- 
plying the absence of a solar-like dynamo in these stars. 
Indeed, no spotted stars earlier than spectral type F8 
(Toff ~6300) have been observed. An exception is the 
class of chemically peculiar Ap and Bp stars, which re- 
tain magnetic fields left over as fossils from the formation 
process or produced by a turbulent dynamo in the core 
(Moss 1972). Objects in this class also exhibit large- 
scale chemical inhomogeneities on their surfaces, which 
are known to induce variability in the light curves of some 
Ap and Bp stars. 

Although observational evidence for spots on high- 
mass stars remains weak, periodic variability in B-type 
main sequence stars attributed to rotational modula- 
tion has been detecte d in high-precision photometry from 
the C oRoT satellite (jPegroote et al.l 120111 iPapics et al. 
2011") as well as with ground-based data (jBriguet et aL 
2001a, b, 2004). The periods of these B stars are be- 
tween 0.3 and 2 days, while their amplitudes are 1-2% 
for the ground-based objects and a much lower ~0.1-1 
mmag for the objects observed with the high-precision 
CoRoT instrument. Perhaps more relevant to the ques- 
tion of HD 31305 are time series photometry from the Ke- 
pler mission indicating that up to 8% of main sequence 
A stars display spot activity (Balona 20111). Notably, 
though, all of these detections were at an amplitude level 
well below 0.5 mmag (see their Fig. 3)- a factor of 30 
lower than the variability we detect in HD 31305. The 
possibility of spots on hot, non-peculiar stars is never- 
theless supported by recent detections of weak magnetic 
fields in a handful of main sequence A stars including 
Veg a (iLignieres et al.l [2001 iPetit et al.l 120101 ) and Sir- 
ius (iPetit et al.l 20111) , as well as several hotter stars 
(|Alecian et al.l 120 lit ). While the paradigm of limited 
magnetic activity in these objects may be weakening, the 
variability characteristics associated with them still does 
not appear to be a good match for the behavior we have 
observed in HD 31305. 

In contrast, the younger Herbig Ae/Be stars are en- 
veloped in a complex and dynamic environment includ- 
ing circumstellar disks, inflowing material, and possible 
outflows, which might provide an explanation for the 



aperiodic components of HD 31305's light curve. These 
objects may also possess short-lived magnetic fields as- 
sociated with the protostellar collapse and accretion 
processes. Such fields have been invoked to explain 
prominent X- ray emission observed in Herbig Ae stars 
(jSwartz et al.l[2005. ') and could be strong enough to gen- 
erate surface spots responsible for the observed periodic 
variability. However, an analysis of the light curve s of 
230 Herbig Ae/Be stars ()Herbst fc ShevchenkollT999l ) re- 
vealed that while nearly all objects were variable, few 
displayed periodicities, on timescales up to 30 days, with 
amplitudes greater than several percent. An additional 
possibility involving magnetic field generated variability 
is that we are observing photometrically not the rota- 
tion of the stellar surface, but the passage of dust clouds 
in the inner disk. The location of material orbiting at 
a Keplerian rotation period of 3.0 days around a young 
AO star with an assumed mass of 2.4 Mq is ~0.05 AU. 
This is an order of magnitude closer to the star than 
the location implied by the 35 -390 K dust temperature 
estimates of lOudmaiier et all ()1992|) and Mooley et al. 
(2013, in preparation). Therefore, for the variability to 
be associated with a disk process, we speculate that co- 
pious material would have to be spiraling in toward the 
host star. Without further information about the nature 
of HD 31305's inclination and infrared excess, it is diffi- 
cult to distinguish this scenario from periodic variability 
generated at the stellar surface. 

A final option is that HD 31305 itself is not variable, 
but rather has a spotted lower mass companion con- 
tributing to the overall flux vari ations and producin g the 
X-ray emission, as proposed bv lArzner et al.l (|2007[ ). We 
checked the Keck Observatory Archive and discovered 
an exposure of this star taken with the NIRC2 adaptive 
optics imaging system as part of engineering tests. The 
Kp image reveals that there is indeed a companion lying 
0.52" away, at a position angle of ~128.5°. We evaluate 
the possibility that this object is the source of variability 
in our MOST observations by considering its brightness 
and inferring the implied variability amplitude. Based 
on the NIRC2 image we estimate the companion to be 
1.3 mag fainter than the primary AO star at the Kp band. 
Since a 0.5" separation indicates a significant probabil- 
ity of association, we consider the case that the two stars 
are bound members of th e ~3 Myr Taurus- Auriga a sso- 
ciation. The TA-DA tool ()Da Rio &: Robbertol[20H en- 
ables us to compute isochrones, infer a mass or tempera- 
ture for the companion, and calculate its expected bright- 
ness in the MOST band, assuming no extinction. In- 
putting AKp=1.3, we find an expected companion tem- 
perature of '^4900 K, or K2 spectral type. Integrating 
a K2 model spectrum over the MOST band from 350- 
750 nm, we predict the companion to be 3.7 mag fainter 
at optical wavelengths. Because of MOSTs 3" pixel size, 
the light curve of HD 31305 includes all of the flux from 
this neighboring star. Thus the lower mass star con- 
tributes just over 3% of the total flux. To cause the 
observed ~1.5% variations, the companion's brightness 
would have to vary by 47%, or 94% peak to peak. This 
value is at the upper end (i.e., <1% level) of the stel- 
lar activity amplitude dis tribution for simi lar tempera- 
ture stars (see Fig. 3 of iBasri et all I2010D , and is ex- 
ceedingly rare among K-type periodic pre main sequence 
stars (e.g.. iStassun et al... 1999, .Cohen et al...2004 ). Fur- 
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thermore, it is difficult to conceive of a spot large enough 
to cause such a substantial variability amplitude. Thus a 
spotted companion does not provide a likely explanation 
for the light curve behavior in our MOST time series. 

We conclude that the variability may come from the 
AO star itself. Given the lack of previous detections of 
periodic variability in main-sequence A stars, we believe 
the most likely scenario for HD 31305 is that it is rela- 
tively young and displays a heretofore undetected type of 
variability, which we tentatively associate with its dusty 
excess or a magnetic field. 

5.3. The origin of periodic variability in SU Aurigae 

Our MOST observations SU Aur have resulted in the 
first detection of clear periodic behavior in this star, at 
a period of 2.66 days. Although SU Aur has been mon- 
itored for years, analysis of its light curve has focused 
on the prominent fading events that led to its classifi- 
cation as an UXOR object. There have been a number 
of periodicity claims (see Table 2), but most were based 
on spectroscopic or X-ray data, and many of the periods 
were noted to have large uncertainties or marginal signif- 
icance. We attempted to recover the 2.66 d periodicity 
in publicly available data from the T Tauri photometry 
database of lHerbst et al.l (|1994[ ) but failed to detect any 
peaks at the corresponding frequency (0.38 d~^) in the 
periodogram. It is conceivable that the combination of 
sparse sampling and systematic dimming events in that 
dataset may have masked transitory periodic behavior. 
Higher cadence monitoring is necessary to determine the 
fraction of time for which periodic variability manifests 
in SU Aur's light curve. 

As pointed out in Section 4.2, the vsini value of SU Aur 
is ~66.2 km s~^ , with an uncertainty of ^^4.6 km s~^ 
([Hartmann et al.l 11986*). If th e stellar inclination is in 
line with that of the disk (62°: IAkeson et al.l[200l . then 
the expected equatorial velocity is 75±7 km s~^. In this 
case, the only way to infer an rotation period as long as 
the one we have measured (2.66 d) is if the stellar radius 
is at least ^3.6 Rq. If we discount the disk inclination 
and raise i to 90°, then the radius may be as low as 
3.2 Rq. 

These v alues stand in contrast to the i? = 2.75±0.25 
derived bvl DeWarf et al . (2003) from sp e ctral e nergy dis- 
tribution (SED) fitting. iCohen et al.l (|1989f ) reported 
a much larger radius (3.6 Rq), which we trace to the 
larger value of the bolometri c luminosity derived b y them 
(~13 Lq vs. -6 Lq from IDeWarf et al.l (120031) '). The 
stellar luminosity of SU Aur is difficult to determine 
accurately because of the disk contribution to its flux 
at near-infrared and infrared wavelengths. We have de- 
rived our own luminosity values using photometry listed 
in the SIMBAD database, along with reddening cor- 
rections from a fit to the SED (spectral type 02) and 
band-dependent bolometric corrections. Depending on 
the bands used, the stellar luminosity could lie anywhere 
between 7 and 30 i© . We favor the lower luminosities of 
7-9 Lq, as these result from SED fits to the U, B, and 
V bands only and exclude contaminating emission from 
the disk at longer wavelengths. Taking these systematics 
into account, we estimate the radius of SU Aur to be be- 
tween 2.9 and 3.2 Rq. While we cannot rule out a higher 
luminosity and hence larger radius, the above estimates 



of spectroscopic rotation velocity and the 2.66 d period 
suggest that these values are unphysical. 

The inconsistency between our measured periodicity 
and the inferred rotation rate of SU Aur leads us to be- 
lieve that we may be witnessing the motion of a dust 
cloud or hot spot connected with this star's inner disk. 
On the presumption that we are not observing starspots, 
but rather observing material in Keplerian orbit about 
the star, we can derive its location by adopting the mass 
of 2.0±0.1 Mq derived bv iDeWarf etall (|200l . We find 
an orbital distance of 7.1 x 10^^ cm, or ~0.05 AU. Intrigu- 
ingly, this is very close to the val ue of the inner disk edge 
reported bv lAkeson et al.l (|2002l ) based on interferomet- 
ric observations of SU Aur (0.05-0.08 AU). We therefore 
conclude that the source of periodic variability observed 
during our MOST observations could be a structural or 
thermal feature on the inner disk edge, or possibly a dis- 
crete cloud of material or biting along magnetic field lines 
(e.g.. lUltchin et aLlll997| ) just interior to the disk. 

Measurement of the rotation period of the inner disk 
also has implications for the angular momentum evolu- 
tion of young stars, many of which are believed to be 
mag netically locked to their disks for at least a few Myr 
(see lCaulev et al.ll2012l for a recent summary). Our in- 
ference of a star rotating significantly faster (P ~1.86 
days, based on the vsim and i values adopted above) than 
the inner disk confronts this idea. The magnetic field 
lines connecting these two regions would likely be highly 
distorted by this velocity shear, and the properties of ac- 
cretion flows could be quite different fro m those predicted 
by m odels involving corotating disks (jRomanova et al.l 
120081 ). Alternatively, we could be observing disk mate- 
rial that is located farther out than the inner edge, and 
thus beyond the corotation radius. 

5.4. Implications for future time series observations of 
young stars 

We have seen that full characterization of young star 
variability requires a combination of high photometric 
precision, short cadence, and sufficiently long time base- 
line that is challenging to arrange. The appearance of 
both periodic and aperiodic variability suggests that it 
may be difficult to infer rotation periods with sparsely 
sampled data, as is usually the case from the ground. 
The large range of rotation periods reported for some of 
our targets underlines this problem. 

In the cases of SU Aur and HD 31305, we have fur- 
ther suggested that periodic variability may not repre- 
sent the stellar rotation period at all, if it is instead tied 
to the surrounding circumstellar disk. We therefore high- 
light a need for more comparisons of usini measurements 
with photometric periods detected in young stars with in- 
frared excesses. We suggest that these measurements be 
carefully reviewed for stars displaying infrared excesses 
or spectroscopic accretion signatures suggestive of disks. 
Rotation period distributions for these objects typically 
have been derived from ground-base d data, and in som e 
clusters display bimodal structure (jHerbst et al.l [200^ . 
Periodicities associated with regions many stellar radii 
above the star's surface may contaminate rotation rate 
samples and thus bias our view of angular momentum 
evolution if not all young cluster stars are locked to their 
disks. 

Different observational complications arise for aperiod- 
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ically variable targets. Our timescale analysis has pro- 
vided guidelines for the time sampling rate in photomet- 
ric future campaigns. For objects with aperiodic light 
curve behavior, our results suggest that data need not be 
taken more frequently than every hour to characterize the 
stochastic components of variability in young stars at 1% 
precision. For higher signal-to-noise observations such as 
those presented here, cadences as short as five minutes 
can be used to probe lower amplitudes of the flicker noise. 
In addition, we find a wide range of coherence timescales 
characterizing the aperiodic variability, from 0.2 to 6.2 
days. However, our very small sample may not be repre- 
sentative of all young stars, especially those in the very 
low mass range below 0.5 Mq. In particular, accret- 
ing brown dwarfs may have shorter characteristic vari- 



ability timescales, as meas ured by autocorrelation (e.g., 
ICodv fc Hillenbrandl[20Tll) . We thus encourage further 
high-precision, well sampled long baseline time series ob- 
servations of young cluster members to expand the size 
and diversity of the variability dataset. Clearly, multi- 
wavelength data as well as high-resolution spectroscopy 
will provide much-needed feedback for model develop- 
ment. 
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